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The Truth of polarization densitites 

So far, we have performed the analysis in a theoretical framework that has been exclusively formu- 
lated in terms of phenomenological models, such as the anharmonic oscillator and the phenomeno- 
logically introduced polarization response function of the medium. In the real world application 
of nonlinear optics, however, we should not restrict the theory just to phenomenological models, 
but rather take advantage over the full quantum-mechanical framework of analysis of interaction 
between light and matter. 



X 




Figure 1. Schematic figure of the ensemble in the "small volume". 

In a small volume V (smaller than the wavelength of the light, to ensure that the natural spatial 
variation of the light is not taken into account, but large enough in order to contain a sufficcient 
number of molecules in order to ignore the quantum-mechanical fluctuations of the dipole moment 
density), we consider the applied electric field to be homogeneous, and the electric polarization 
density of the medium is then given as the expectation value of the electric dipole operator of the 
ensemble of molecules divided by the volume, as 

P^{r,t) = {Q^)/V, 

where the electric dipole operator of the ensemble contained in V can be written as a sum over all 
electrons and nuclei as 

Q = -e ^ f-j + e ^ ZkTk ■ 

electrons nuclei 

The expectation value {Q^) can in principle be calculated directly from the compound, time- 
dependent wave function of the ensemble of molecules in the small volume, considering any kind 
of interaction between the molecules, which may be of an arbitrary composition. However, we will 
here describe the interactions that take place in terms of the quantum mechanical density operator 
of the ensemble, in which case the expectation value is calculated from the quantum mechanical 
trace as 

P^{r,t) = TrmQ^^]/V. 
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Outline 

Previously, in lecture one, we applied the mathematical tool of perturbation analysis to a classical 
mechanical model of the dipole moment. This analysis will now essentially be repeated, but now 
we will instead consider a perturbation series for the quantum mechanical density operator, with 
the series being of the form 

p{t)= ^ + pi{t) + P2{t) +...+ Pn{t) +... 

-^[E(t)Y> -[B(t)]^ -[BW]" 

As this perturbation scries is inserted into the expression for the electric polarization density, we 
will obtain a resulting series for the polarization density as 

oo n 
m=0^ ~X ' m=0 



Quantum mechanics 

We consider an ensemble of molecules, where each molecule may be different from the other 
molecules of the ensemble, as well as being affected by some mutual interaction between the other 
members of the ensemble. The Hamiltonian for this ensemble is generally taken as 

H = Ho + Hi{t), 

where Hq is the Hamiltonian at thermal equilibrium, with no external forces present, and Hi{t) 
is the interaction Hamiltonial (in the Schrodinger picture), which for electric dipolar interactions 
take the form: 

Hi{t) = -Q • E(r,i) = -Qc,Ec,{r,t), 

where Q is the electric dipole operator of the ensemble of molecules contained in the small volume V 
(sec Fig. 1). This expression may be compared with the all-classical electrostatic energy of an 
electric dipole moment in a electric field, V = — p • E(r, t). 

In order to provide a proper description of the interaction between light and matter at molecular 

level, we must be means of some quantum mechanical description evaluate all properties of the 
molecule, such as electric dipole moment, magnetic dipole moment, etc., by means of quantum 
mechanical expectation values. 

The description that we here will apply is by means of the density operator formalism, with 
the density operator defined in terms of orthonormal set of wave functions |a) of the system as 

P = YPa\a){a\ = p{t), 

a 

where Pa ^^"^ the normahzed probabihties of the system to be in state |a), with 

a 

From the density operator, the expectation value of any arbitrary quantum mechanical operator 
O of the ensemble is obtained from the quantum mechanical trace as 

{6)=Tr{pd) = Y,{k\pO\k). 

k 

The equation of motion for the density operator is given in terms of the Hamiltonian as 

= [Ho,p] + [Hi{t),P\ 
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In this context, the terminology of "equation of motion" can be pictured as 

J A change of the density ^ J A change of density 1 J Change of a general 1 

\ operator p{t) in time J I of states in time J \ property (O) in time J 

Whenever external forces arc absent, that is to say. whenever the applied electromagnetic field is 
zero, the equation of motion for the density operator takes the form 

with the solution^ 

p{t) = po = r]exp{-Ho/kBT) 

being the time-independent density operator at thermal equilibrium, with the normalization con- 
stant T] chosen so that Tr(p) = 1, i. e., 

77 = l/TV[exp(-^oABT)]. 
Perturbation analysis of the density operator 

The task is now o obtain a solution of the equation of motion (1) by means of a perturbation 
series, in similar to the analysis performed for the anharmonic oscillator in the first lecture of this 
course. The perturbation series is, in analogy to the mechanical spring oscillator under infiuence 
of an electromagnetic field, taken as 

p{t)= ^ + + p2{t) +...+ Pn{t) +... 

~[B(t)i° -imv -{E{tw -[BW]" 

The boundary condition of the perturbation series is taken as the initial condition that sometime 
in the past, the external forces has been absent, i. e. 

p(-oo) = po, 

which, since the perturbation series is to be valid for all possible evolutions in time of the externally 
applied electric field, leads to the boundary conditions for each individual term of the perturbation 
series as 

Pji-oo) = 0, i = l,2, ... 

By inserting the perturbation series for the density operator into the equation of motion (1), one 
hence obtains 

i/i^(/3o + Piit) + h{t) + ... + Pn{t) + ...) = [Ho, Po + Pi{t) + h{t) + ... + Pn{t) + ...] 

+ [Hi{t), Po + hit) + hit) + ■■■ + Pnit) + •••], 



^ For any macroscopic system, the probability that the system is in a particular energy eigenstate 
V'n, with associated energy E„, is given by the familiar Boltzmann distribution 

Pn = r7exp(-E„/A:Br), 

where r? is a normalization constant chosen so that J2nPn = !> ''^b is the Boltzmann constant, 
and T the absolute temperature. This probability distribution is in this course to be considered 
as being an axiomatic fact, and the origin of this probability distribution can readily be obtained 
from textbooks on thermodynamics or statistical mechanics. 
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and by equating terms with equal power dependence of the appUed electric field in the right and 
left hand sides, one obtains the system of equations 



dt 
dt 
dt 



in^^ = [Ho,Pi{t)] + [Hi{t),po], 



iH^ = [Ho,p2{t)] + [Mt),Pi{t)], 



for the variuos order terms of the perturbation series. In Eq. (2) , we may immediately notice that 

the first equation simply is the identity stating the thermal equilibrimn condition for the zeroth 
order term po, while all other terms may be obtained by consecutively solve the equations of order 
j = 1, 2, . . . , n, in that order. 

The interaction picture 

We will now turn our attention to the problem of actually solving the obtained system of equations 
for the terms of the perturbation series for the density operator. In a classical picture, the obtained 
equations are all of the form similar to 

f^=f{t)p + g{t), (3) 

for known functions f{t) and g{t). To solve these equations, we generally look for an integrating 
factor I{t) satisfying 

7(t)^-7(t)/(t)p=|[7(i)p]. (4) 

By carrying out the differentiation in the right hand side of the equation, we find that the inte- 
grating factor should satisfy 

'-f-nm,, 

which is solved by^ 

/(i)=/(0)exp[- / /(T)rfr]. 

Jo 

The original ordinary differential equation (3) is hence solved by multiplying with the intagrating 
factor I{t) and using the property (4) of the integrating factor, giving the equation 

|[J(t)p]=J(i)5(i), 
from which we hence obtain the solution for p{t) as 

Prom this preliminary discussion we may anticipate that equations of motion for the various order 
perturbation terms of the density operator can be solved in a similar manner, using integrating 



^ Butcher and Cotter have in their classical description of integrating factors chosen to put 
7(0) = 1. 
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factors. However, it should be kept in mind that we here are dcahng with operators and not 
classical quantities, and since we do not know if the integrating factor is to be multiplied from left 
or right. 

In order not to loose any generality, wc may look for a set of two integrating factors Vo(f) 
and Uo{t), in operator sense, that we left and right multiply the unknown terms of the nth order 
equation by, and we require these operators to have the effective impact 

Vo{t){ih^^ - [Ho,Pn{t)]^Uo{t) = ifif^[Voit)Pnmo{t)]- (5) 

By carrying out the differentiation in the right-hand side, expanding the commutator in the left 
hand side, and rearranging terms, one then obtains the equation 



ihjVoit) + Voit)Ho jpnitWoit) + Vo{t)Mt)\it^jUo{t) - HoUo{t) | = 

for the operators Vo{t) and Uo{t)- This equation clearly is satisfied if both of the braced expressions 
simultaneously are zero for all times, in other words, if the so-called time-development operators 
Vo{t) and Uo{t) are chosen to satisfy 

ih^ + Vo{t)Ho = 0, 
ih^-HoUoit) = 0, 

with solutions 

Uoit) = expi-iHot/h), 

Vo{t) = exp{iHot/h) = Uoi-t). 

In these expressions, the exponentials are to be regarded as being defined by their series expansion. 

In particular, each term of the series expansion contains an operator part being a power of the 
thermal equilibrium Hamiltonian Hq, which commute with any of the other powers. We may easily 
verify that the obtained solutions, in a strict operator sense, satisfy the relations 

Uoit)Uo{t') = Uo{t + t'), 

with, in particular, the corollary 

Uo{t)Uoi-t) = Uo{0) = 1. 

Let us now again turn our attention to the original equation of motion that was the starting point 
for this discussion. By multiplying the nth order subequation of Eq. (2) with Uo{—t) from the left, 
and multiplying with Uo{t) from the right, we by using the relation (5) obtain 

ifij^l^Uo{-t)Pnmo{t)^ = M-t)[Hl{t),Pn-l{t)\Uo{t), 
which is integrated to yield the solution 

M-t)Pn{t)Mt) = ^f Uo{-r)[Hi{r),Pn-i{r)]Uo{r)dT, 

where the lower limit of integration was fixed in accordance with the initial condition cx)) = 0, 

n = 1, 2, In some sense, we may consider the obtained solution as being the end point of this 

discussion; however, we may simplify the expression somewhat by making a few notes on the 
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properties of the time development operators. By expanding the right hand side of the solution, 
and inserting tJo(T)ijQ{—T) = 1 between Hiij) and p„_i(r) in the two terms, we obtain 

C/o(-t)p„(i)C>o(i) = ^ /* M-r)[Hi{r)Pn-i{T) - Pr,_iiT)H,{T)]Uo{T) dr 

in y_oo ^ , ' 

=1 

- 4 /* Uoi-T)pn-i{T) Uo{r)Uoi-r) Hr)Mr) dr 



f [Uo{-T)Hi{T)Uo{T)M-r)Pn-i{r)Uo{r)]dT, 

'■n J-oo '• V ' ^ V/ ' 



and hence, by introducing the primed notation in the interaction picture for the quantum mechan- 
ical operators, 

pUt) = lfo{-t)Pnmo{t), 

H[{t) - Uo{-t)Hi{t)UQ{t), 

the solutions of the system of equations for the terms of the perturbation series for the density 
operator in the interaction picture take the simplified form 



P'n{t) = l^f JH[{T),p'n-,{r)]dT, 



n=l,2,..., 



with the variuos order solutions expressed in the original Schrodinger picture by means of the 
inverse transformation 

Pn{t) = Uo{t)p'^{t)Uo{-t). 



The first order polarization density 

With the quantum mechanical perturbative description of the interaction between light and matter 
in fresh mind, we are now in the position of formulating the polarization density of the medium 
from a quantum mechanical description. A minor note should though be made regarding the 
Hamiltonian, which now is expressed in the interaction picture, and hence the electric dipolar 
operator (since the electric field here is considered to be a macroscopic, classical quantity) is given 
in the interaction picture as well, 

Hiir) = Uo{-t) l-QaEair)] Uo{t) 

^ V ' 

=Hi(r) 

= -Uo{-t)Q^Uo{t)E„{t) 
where Qa (t) denotes the electric dipolar operator of the ensemble, taken in the interaction picture. 
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By inserting the expression for the first order term of the perturbation series for the density 
operator into the quantum mechanical trace of the first order electric polarization density of the 
medium, one obtains 



pW(r,t) = iTr[pi(i)Q^] 



V 

1 - 



= -Tr 
V 



{Uo{t) ^ J [Hi{r),Po] drUoi-t)) Q, 



=pi(«) 



=Pi(t) 

Eij,{t) is a classical field (omit space dependence r), 
[H{{T),h] = [-QaiT)E^{T),po] = --E„(r)[Q„(r),/5o] 

TrhJrJi'] ( 

Vih 

= {Pull out Ea{T) and the integral outside the trace} 

Vih 



(r ) Tr{C/o (i) [Q„ (r ) , po] Uo{-t)Q^} dr 
= I Express ^^^(t) in frequency domain, i?c(r) = J Ea{u)) exjp{—iu)T) dio^ 

POO pt 

^~VihJ J Ec.{uj)Ti{Uo{t)[Qa{T),po]Ua{-t)Q^,}exp{-iwT)dTduj 



-oo «/ — oo 

= {Use exp(— zojt) = exp(— ia;t) exp[— ^ct;(T — t)]} 

(w) Tv{Uo (i) [Qa (r) , M M-t)Qf,] 



1 



CO J —CO 



X exp[— zu;(t — f)] dr exp(— iwt) 



£0 / xLa(-'^;^)-^a(w)exp(-iwt)dw, 



OO 



where the first order (linear) electric susceptibility is defined as 



xL^2(-w;w) 



1 . . 

= TTTz / Tr{i7o(t)[0„(r),po]C/o(-i) Q^} expHa;(r - t)] dr 

= [Qa.(r-l),M 

{Expand the commutator and insert UQ{—t)Uo{t) = 1 ~) 
in middle of each of the terms, using [Uoit),po] = f 

= ^ / Tr{[Q„{T-t),po]Q^.}eM-Mr-t)]dr 

= < Change variable of integration t' = t — t; ■ ■ ■ dr ^ / ••• dr' > 

I J —oo J —oo J 

= T7TZ / Tr{[g„(T'),/3o]Q^}exp(-iu;T')dr' 

=Tr{po[(3^,Q=(T')]} 

= {Cyclic permutation of the arguments in the trace} 

1 - ^ 

= TTTT Tr{po[Q/:i,Q«(-r)]}exp(-iwr)(iT. 
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